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Abstract 

We investigate interactions between an elastic membrane and a substrate charac- 
terized by quenched positional disorder in the height function. We show that the posi- 
tional disorder transforms the standard secondary DLVO minimum into two separate 
states: the hovering state characterized by a planar membrane at a finite separation 
from the interface and a pinned state where the membrane follows closely the asperi- 
ties of the substrate and is as a consequence quite corrugated. The transition between 
the two states is continuous and depends on the parameters of the underlying DLVO 
potential as well as the parameters describing the quenched height-height correlation 
function of the substrate. 

pacs87.16.DgMembranes, bilayers, and vesicles pacs68.15.-|-eLiquid thin films 

In the DLVO theory one usually assumes that when a membrane or in general an elas- 
tic manifold interacts with a rigid substrate, the membrane is modelled as flat and the 
substrate is envisioned as featureless This leads in a straightforward manner to the 
secondary minimum of the DLVO theory and to an equilibrium spacing between the mem- 
brane and the substrate. Taking into account the elastic degrees of freedom of a membrane 
can alter this picture drastically, leading to the emergence of Helfrich interaction which 
may induce an unbinding transition of the membrane. This transition has been clearly 
shown to be a consequence of the interplay between elastic fluctuations of the membrane 
and DLVO interactions between the membrane and the substrate 0. In the present 
work we are motivated by recent experiments on deposited lipid multilayers on atomically 
smooth vs. rough surfaces |^. 

These experiments reveal differences in equilibrium lamellar spacings of lipid multilayers 
in proximity of a substrate that seem to correlate with molecular roughness of the substrate 
H . In order to lay ground for understanding effects of this type it is desirable to relax the 
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model constraint of a featureless substrate and treat it supposedly more realistically as 
exhibiting (quenched) disorder in the height function. We will analyze the consequences 
of this new model of a substrate interacting with a membrane and show that it modifies 
the simplified DLVO conclusions in the sense that the original secondary DLVO minimum 
now splits into two separate states, characterised by the way the membrane is (de) coupled 
to the substrate. We call these states the hovering state and the pinned state. The former 
one is characterized by a membrane in a flat configuration hovering a certain finite distance 
above the substrate and is directly related to the standard DLVO secondary minimum. The 
latter one is characterized by a corrugated membrane that basically follows the asperities 
of the quenched positional disorder of the substrate and is altogether missing from the 
standard DLVO theory. We stress that our description is unifying in the sense that it 



combines aspects of mean-Held theory 
and a straight-forward generalization 14 



|6|, which focuses on the nature of the pinned state; 
of Li and Kardar's |^] Gaussian fluctuation theory, 
which allows for a hovering state determined by the parameters of the DLVO interaction 
potential. Specifically, we find that for sufficiently strong disorder, the membrane indeed 
prefers the mean-field pinned state considered by Swain and Andelman. Because of the 
increased role of substrate disorder fluctuations, in weakly disordered systems the properly 
self averaged free energy can give way to a hovering state where the membrane is depinned, 
residing in a minimum related to but different from the usual DLVO minimum, as predicted 
by Gaussian fluctuation theory. 

We introduce our approach by first specifying the model of the membrane and the 
substrate that interact via a DLVO potential with (for matters of convenience) hydration 
and van der Waals terms. Formally our analysis owes a lot to investigations of effects of 
the disorder on polymerized membranes that have been analyzed in a variety of contexts 
p. It is based on the assumption of a Gaussian substrate height distribution function 
and the application of the 1/d expansion method Q to evaluate the complicated partition 
function stemming from the replicated Hamiltonian. 

The Hamiltonian of a flexible membrane described in the usual Monge parameterization 
u{p), where p = {x,y) is a 2D coordinate vector, above a substrate characterized by the 
height function C(p)) where C(p) is a quenched disorder field describing the profile of the 
substrate, is composed of three contributions. First of all we have the elastic energy of the 
membrane, then the interaction free energy per unit surface area between the membrane 
and the substrate V{u{p) — C(p)) of the general DLVO form, and finally the free energy 
contribution of an external force per unit surface area vr pushing the membrane towards 
the substrate. The total free energy or equivalently the mesoscopic Hamiltonian thus 
assumes the form 

nHp)] = h IJd'pd'p' JC{p, p'MpHp') + Jd'p V{u{p) - C(p)) + Jd'p n{u{p) - C(p)). 

Usually one takes for the elastic part the classical Canham - Helfrich - Evans ansatz 
}C{p, p') = ^'^KcS'^{p — p'), where Kc is the bending modulus of the membrane, while the 

^Strictly speaking what we use is actually a l/{d — 2) expansion. 
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interaction part is composed of the attractive and repulsive DLVO ingredients. Since the 
interaction part of the Hamiltonian is in general non-linear it is convenient to introduce 
the following new variable B{p) = {u{p) — C{p))'^ at every p via a functional constraint |^] 

6n[B{p),g{p)] = '^Jd'pg{p){{n{p)-ap)r-B{p)). 

Clearly the auxiliary field g{p) plays a role akin to the self-energy part of the Green 
function. The two auxiliary fields B{p) and g{p) just introduced, play the role of fixing 
the local constraint (n(p) — C(p))^ = ^{p)- They also show up in the partition function 
where one has to eventually take the trace over these auxiliary fields together with u{p). 
The partition function can thus be written as 

Z[C(p)] = j Vu{p)VB{p)Vg{p) e'mHP)]+mB{p),9m) . 

The average over the quenched disorder distribution, assuming the self-averaging property 
of the free energy, is defined to be of the form 



{■■■) = ] vap){. . .)v[ap)i 

where by assumption the disorder probability distribution function V{({p)) that charac- 
terises the quenched disorder in the height function of the substrate is given by a Gaussian 
ansatz 

v[ap)] = exp{-^jjd'pd^p' g{p,p')ap)ap')). 

The free energy, after being averaged also over the quenched disorder distribution, is 
obtained finally as 



z-[C{p)] - 1 



JF= -kT log Z[C(p)] = -kT lim ^-^^ . (1) 

n — >o n 

With these preliminaries the free energy Eq. || can be evaluated via the standard Edwards- 
Anderson replica trick Q , where the replicated Hamiltonian TCn is composed of the repli- 
cated elastic term 

hU=oIId'pd'p' lC{p,p')u,{p)u,{p'), 

where i is the index of the replica, the replicated constraint on the variable Bi{p) that 
now reads 

I E-=o fd'p 9^{P) {{ui{p) - Cip)? - B,{p)) , 
and finally of the replicated interaction and external "source" terms 

n „ n „ 

d'p V{B,ip)) + d^P 7Tiu,{p) - C(P)). 

1=0 1=0 
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Since we now have Gaussian integrals over the variables Ui{p), C(p) we can evaluate them 
explicitely, while the functional integrals over the auxiliary fields Bi{p),gi{p) can be evalu- 
ated on the saddle-point level, with the proviso that there is no replica symmetry breaking. 
This constitutes the essence of the 1/d expansion. 

After performing all the indicated integrations and taking the n — > limit, the free 
energy can be obtained as a sum of the mean-field part and a fluctuation part. The 
mean-field part is 

To = J(fpiTUoip) + ^JJ(fp(fp'xiP,p')uo{p)uoip') 

where the mean-field uq(p) is obtained via minimization of J-q with 

x{p,p') = }Cip,p') + 5\p-p')g{p). 
The fluctuation part of the free energy is concurrently obtained as 

T = ^ Jd^pg-\p,p)g{p)-^ JJd^pd'p'x-Hp,p')g-Hp',p)g{p)g{p')- 

- \!d'pg{p)B{p) + fd'pV{B{p)). (2) 

On the \/d expansion level the auxiliary fields contribute only at the saddle point. The 
saddle point equations are now obtained simply by minimizing Eq. ^ with respect to B{p) 
and g{p). We will not reproduce the general rather awkward form of these equations 
but will concentrate on a rather particular solution characterized by uq = const, and 
B = const. ^ implying also g = const.. In addition to this we will limit ourselves to the 
conceptually most interesting case of vanishing external confining force, i.e. vr = 0. 

Assuming that the system is homogeneous in the p plane we can introduce the Fourier 
transforms of all the relevant quantities that allow us to write the mean-field equations in 
a rather simple form 

uox{Cl = 0) = no5 = 0. (3) 
The saddle point for B is obtained straightforwardly as 

dV{B) _ 1 

dB 

while the saddle point for g can be derived in the form 

\2 



B = ul + kTY' 



Q 



g(Q)(/C(Q)+5)2' 



This set of equations has two fundamentally different solutions describing the state of the 
elastic membrane interacting with a disordered substrate. 

^Not to be confused with the mean-field approach introduced in [^. The term mean-field is used here 
strictly as it pertains to the 1/d expansion. 
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1. a hovering state, with uq ^ 0, g = 0, characterized by 
dV{B) 



=0 with B = ul + kT^g-\Q). (4) 

Q 

2. and a pinned state, with uq = 0, g 0, characterized by 

^y^-^a with B-kT^ ^(Q)^ 

Obviously, in order to progress we have to assume a certain form for the quenched disorder 
correlation function. As a simplest approximation we take a spatially short range coupling 
Q{p,p') = GS'^{p — p'), implying ^(Q) = G. Thus we have in the hovering state 

^^^ = with ul = B-Bc, (6) 

where we introduced B^ = fer^Q^-i(Q) = that obviously depends on the 

upper wavevector cutoff in the Fourier space. The first of the above equations determines B 
as a function of the parameters of the DLVO potential {e.g. Hamaker constant, hydration 
interaction strength etc.). The second one gives the dependence of uq on these parameters. 

In order to understand the physical nature of the two phases, we evaluate the average 
of the separation between the substrate and the membrane that can be obtained as 



Jd^p<{u{p)-ap))> = uo. 



Clearly uq quantifies the disorder averaged separation between the membrane and the 
substrate. For a finite uq the membrane hovers a finite separation away from the substrate. 
In order to characterize the hovering state further we evaluate the average square of the 
deviation from a unit normal to the membrane in the z direction given to the lowest order 
as 5n(p) Vj_ • u{p). 

Thus in the hovering state with ^ = not only is the membrane decoupled from the 

substrate and hovers above it, but is also flat on the average since the mean disorder 
averaged squared deviation from the normal, p^, equals zero at a finite separation uq (see 
figure). Obviously the hovering solution exists only for B > B^. The hovering line in the 
"phase diagram" thus ends at the value of the interaction parameters where the solution 
of ^g^^ = also satisfies B = Be- uq thus behaves as an order parameter of a second 
order phase transition, and B behaves as the temperature. 
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Furthermore the surface density of the adhesion free energy ([/) in the hovering state is 
jiven by 



C/ = mm = V{B) = V{ul + B,). 



This obviously differs from the standard secondary minimum of the DLVO theory, deter- 
mined from = 0, with U = V{uq). The free energy in the hovering state is thus not 
given by the value of the DLVO potential at the minimum V{uq)^ but at a (much smaller) 
value of V{uq + B c), since is a decreasing function of its argument above its minimum. 

In the hovering state the membrane is obviously in close proximity ("contact") of the 
substrate only for a fraction of its total surface area. This effective contact area [Sc) can 
be estimated from 

S, _V{ul + B,) 



S y(n^) 

In view of the discussion presented above, this ratio is smaller then one if ^ is a decreasing 
function of its argument above its minimum value. 

For the pinned state, with /C(Q) = Q'^Kc, we obtain the following expression to the 
lowest order in the wavevector cutoff Qmax 



- 29 with ^-^c-T6GV^' 



where Be is defined in the same way as in the hovering state. The adhesion free energy in 
the pinned state, defined in complete analogy with the hovering state, is given by 

In the pinned state the membrane is thus coupled to the substrate and follows it closely, 
being always in its close proximity since = (see figure). It thus exhibits a very 
corrugated configuration. This follows again from Eq. |^ since in the pinned state ^ 0. 
We note that for non-zero external driving force tt the system is always in the pinned 
state. 

The adhesion energy in both states is nowhere in general equal to its DLVO counterpart. 
Both, the hovering as well as the pinned states carry in the free energy the signature of 
the substrate disorder in the height-height correlations. Only in the limit of vanishing 
disorder, or in the language of our model as G — > oo, does the hovering state approach 
the DLVO secondary minimum in a continuous fashion, while the pinned state simply 
disappears. 

In order to gain further insight into the nature of the hovering and the pinned states 
we investigate the phase diagram for a particular typical choice of the DLVO interaction 
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potential augmented by the Helfrich undulation interaction. For the sake of simplicity 
we assume that the DLVO part is given by the sum of the hydration and van der Waals 
interactions while the Helfrich interaction is assumed to have the same form as between a 
flexible membane and a flat substrate |l|, ^, thus 

where A is the magnitude and ^ the range of the hydration interaction, and treat the 
Hamaker constant H as a variable tuning parameter. The above form of the total in- 
teraction energy between a membrane and a substrate would be strictly valid only for a 
membrane fluctuating near a flat substrate. However Swain and Andelman note Q that 
the form of the Helfrich interaction should not be far off from the one in the above expres- 
sion even for a rough substrate. We assume this is the case when we analyze the phase 
diagram for a membrane near a disordered substrate. This assumption would however 
have to be tested via a more sophisticated and hopefully more accurate approach. 

Without the disorder the interaction potential Eq. |8| leads to an unbinding transtion at 
the critical value of the Hamaker constant equal to H^. = "^-^gg^-^- When H approaches 
this value the secondary DLVO minimum is displaced towards infinity and we have a 
continuous unbinding of the membrane. Adding disorder to this scenario we instead obtain 
a modified phase diagram as presented on the figure. There is now a window corresponding 
to the hovering state (the former DLVO secondary minimum) in between the pinned and 
the unbound state of the membrane. The dimensions of this window depend on the value 
of Q that characterizes the intensity of the substrate disorder. In the limit of no disorder 
Q — > oo with Be — > the hovering window is expanded to the whole axis H > He and 
is tranformed back into the standard DLVO secondary mimimum. 

The disorder usually does not figure in the theories of membrane substrate interactions. 
The present work is basically a plea for a change of this perspective. There obviously 
exist phenomena, where ignoring the disordered nature of the substrate does not, even 
qualitatively, lead to the correct physical picture. In this sense the DLVO theory has to 
be ammended. 
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Figure 1: Upper left: Schematic dependence of uq on the auxihary field B. The hovering 
state exists only for B > Be- The critical value Be depends on the disorder characteristics 
of the substrate and the microscopic wavevector cutoff. The hovering - pinned transition 
is obviously second order in B, which thus plays the role of the temperature. Upper right: 
Schematic phase diagram of a membrane as a function of the Hamaker coefficient; V{u) is 
assumed to be equal to the sum of hydration, van der Waals and Helfrich terms Eq. ^. The 
substrate disorder introduces a hovering (H) state window into the phase diagram, located 
between the pinned (P) and the unbound (U) phases. Two possible (H-U) boundaries are 
indicated by the HI and H2 lines. The size of this window depends on the strength of the 
Hamaker coefficient and the value of G, characterizing the disorder effects. In the case of 
no disorder (B — > 0) the critical value of the Hamaker coefficient He marks the unbinding 
transition of the membrane. Bottom: A schematic representation of the hovering (H) and 
the pinned (P) state of the membrane in proximity of a disordered substrate. The hovering 
state is characterized by a slowly varying average separation between the membrane and 
the substrate that does not follow closely the local corrugations of the substrate. In the 
pinned state however, the membrane follows closely the asperities of the substrate. 
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